Anomalous Periodicity of the Current-Phase Relationship of Grain-Boundary 
Josephson Junctions in High-T c Superconductors 
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The current-phase relation (CPR) for asymmetric 45° Josephson junctions between two d-wave 
superconductors has been predicted to exhibit an anomalous periodicity. We have used the single- 
junction interferometer to investigate the CPR for this kind of junctions in YBa2Cu307-i thin films. 
Half-fluxon periodicity has been experimentally found, providing a novel source of evidence for the 
d-wave symmetry of the pairing state of the cuprates. 



There is growing evidence in favor of the d x 2_ y 2-w&ve 
symmetry of thenpairing state of the high-temperature 
superconductors.u An unconventional pairing state re- 
quires the existence of zeros of the order parameter in 
certain directions in momentum space. Thermodynamic 
and spectroscopic measurements do indeed suggest their 
existence, but by themselves they do not exclude conven- 
tional s-wave pairing with nodesu Direct evidence for 
the d-wave pairing state is provided by phase-sensitive 
experiments, which are based on the Josephson effect.0 
Quite generally, the current-phase relationship (CPR) of 
a Josephson junction, J(f), is an odd periodic function 
of (p with a period 2ir.a Therefore 1(f) can be expanded 
in a Fourier series 



I (ip) = I\ sin (f + 1% sin 2ip ■ 



(1) 



In the tunnel limit we can restrict ourselves to the first 
two terms in Eq. (Q). Since the order parameter is 
bound to the crystal lattice, I((p) of a weak link de- 
pends on the orientation of the <i-wave electrodes with 
respect to their boundary. The existing phase-sensitive 
experiments exploit possible sign changes of Ji between 
different geometries. In this Letter we present a new 
phase-sensitive experimental test of the symmetry of the 
pairing state of the cuprates. Namely, in certain geome- 
tries, the I\ term should vanish by symmetry. In such 
cases, the CPR should exhibit an anomalous periodicity. 

Let us analyze the angular dependence of in a junc- 
tion between two macroscopically tetragonal <i-wave su- 
perconductors. As emphasized in Ref. [|, also heavily 
twinned orthorhombic materials such as YB^CusOy-z 
belong to this class, if the twin boundaries have odd sym- 
metry. We consider first an ideally flat interface between 
the superconducting electrodes. Let B\ (0 2 ) denote the 
angle between the normal to the grain boundary and the 
a axis in the electrode 1 (2), see Fig. 1. If we keep only 
the lowest-order angular harmonics, the symmetry of the 
problem dictates thafn 

h = I c cos 20i cos 26*2 + I s sin 20i sin 20 2 . (2) 



The coefficients I c , I s are functions of the barrier 
strength, temperature T, etc. The I 2 term results from 
higher-order tunneling processes and we neglect its weak 
angular dependence. It is seen from Eq. (|^) that the cri- 
terion for the observation of an anomalous CPR, I\ = 0, 
is realized for an asymmetric 45° junction, i.e. a junction 
with 0\ = 45° and 9 2 = 0. For an interface which is not 
ideally flat, 9i = 9i(x) are functions of the coordinate x 
along the junction. Ii = remains valid even in this 
case, if the average values (9i(x)) = 45° and (#2(2:)) = 0. 




FIG. 1. Sketch of the interferometer (not in scale). 

The I 2 term is present also in weak links based on 
conventional s-wave superconductors but for all known 
types of weak links l-fe/iil < 1- For instance, for a tunnel 
junction <C 1- For a SNS junction, / oc sin^/2 

at T = 0,Ej and the Fourier expansion Eq. (|l]) leads to 
I2/I1 = —2/5. Therefore, a possible experimental ob- 
servation of l/2/^i I S> 1 in an asymmetric 45° junction 
provides direct evidence of <i-wave symmetry of the pair- 
ing state in the cuprates. 

We have investigated the CPR of YBa2Cu307_ x 
thin film bicrystals with asymmetric 45° [001]-tilt 
grain boundariea3~13 as sketched in Fig. 1, using a 
single-junction interferometer configuration in which the 
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Josephson junction is inserted into a superconducting 
loop with a small inductance L. In a stationary state 
without fluctuations, the phase difference tp across the 
junction is controlled by applying external magnetic flux 
<& e penetrating the loop: <p — (p e — f3f(ip). Here (p e — 
27r$ e /$o is the external flux normalized to the flux quan- 
tum $ = 2.07 x 1CT 15 Tm 2 . The CPR is written as 
I(ip) = Iof(<p), where Iq is the maximal Josephson cur- 
rent. j3 = 2ttLIo/$q is the normalized critical current. 
In order to obtain the CPR for the complete phase range 
— 7r < (p < 7r the condition f3 < 1 has to be fulfilled, 
because for (3 > 1 the curve <p{<p e ) becomes multivalued 
and there are jumps of (p and a hysteresis for a sweep 
of ip e . Following Ref. [| we express the effective induc- 
tance of the interferometer using the derivative /' with 
respect to ip as L; mt — L[l + 1/ f'(ip)]. The inductance can 
be probed by coupling the interferometer to a tank cir- 
cuit with inductance Lt, quality factor Q, and resonance 
frequency w .Ij External flux in the interferometer is pro- 
duced by a current Idc + Irf in the tank coil and can be 
expressed as tp e = 2ir(Id c + I r f)M /$ — ipd c + Lp r f, where 
M 2 = k 2 LLx, and k is a coupling coefficient. Taking into 
account the quasiparticle current in the presence of volt- 
age V across the junction the phase difference is given 
by the relation ip — ipdc + <p>rf — Pf(p) ~ 2irT(tp)V/$o, 
where r(ip) = L/Rj(ip) and Rj(<p) is the resistance of 
the junction . In the small-signal limit ip r f -C 1 and in 
the adiabatic case cot <C 1, keeping only the first-order 
terms, the effective inductance L e ff of the tank curcuit- 
interferometer system reads 



L 



eff 



L T 1 - k 2 



L 



U 



Ln 



1 



k 2 f3f'(<p) 
1 + Pf'(<p) 



Thus the phase angle a between the driving current and 
the tank voltage U at the resonant frequency of the tank 
circuit loq is 



t&na(cp) — 



k 2 Q(3f(y) 
1 + 



(3) 



Using the relation [1 + (3f'((p)]dip — dipdc valid for iprf -C 
1 and lot -C 1, one can find the CPR from Eq. (m by 
numerical integration. 

The advantage of the measurement of the CPR of an 
asymmetric 45° junction with respect to the by-now stan- 
dard phase-sensitive tests of pairing symmetry based on 
the angular dependence of I\ is twofold. First, it avoids 
the complications of the analysis of experiments caused 
by the presence of the term I s a Second, a flux trapped 
in the SQUID does not invalidate the conclusions about 
the ratio IJ2/J1I and hence about the pairing symmetry, 
while this is not the case in standard phase-sepsitive tests 
of the d-wave symmetry of the pairing state.li9 

The films of thickness 100 nm were fabricated us- 
ing standard pulsed laser deposition on (001) oriented 
SrTiC>3 bicrystalline substrates with asymmetric [001] tilt 
misorientation angles 45° ±1°. They were subsequently 
patterned by Ar ion-beam etching into 4x4 mm 2 square 



washer single-junction interferometer structures (Fig. 1). 
The widths of the junctions were 1 -j- 2 fxva. The washer 
square holes had a side-length of 50 fim. This geometry 
of the interferometer gives L ss 80 pH. The resistance of 
the junction is higher than 1 ft and the condition for the 
adiabatic limit lot <gc 1 is satisfied. For measurements of 
a(tpdc), several tank circuits with inductances 0.2 ~ 0.8 
/iH and resonance frequencies 16 -j- 35 MHz have been 
used. The unloaded quality factor of the tank circuits 
70 < Q < 150 has been measured at various tempera- 
tures. The coupling factor k was determined from the 
period Aide of a(Id c ) using MAIjc = $o- Its value var- 
ied between 0.03 and 0.09. The amplitude of I r f was 
set to produce the flux in the interferometer lower than 
0.1$ . 

The measurements have been performed in a gas-flow 
cryostat with a five-layer magnetic shielding in the tem- 
perature range 4.2 < T < 90 K. The experimental setup 
was calibrated by measuring interferometers of the same 
size with 24° and 36° grain boundaries. We have stud- 
ied 5 samples, out of which sample No. 1 exhibited the 
most anomalous behavior. Samples Nos. 2,3 were less 
anomalous and the remaining two samples had high crit- 
ical currents and their I(ip) was conventional. In Fig. 2 
we plot the phase angle a as a function of the dc current 
Idc for samples Nos. 1,2. The data for the 36° junction 
is shown for comparison. Note that at T — 40 K the 
periodicity of a(Idc) is the same for all samples. 
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FIG. 2. Left panel: Phase angle between the driving current 
and the output voltage measured for the sample No. 1 at dif- 
ferent temperatures as a function of the dc current Idc- The 
curves are shifted along the y axis and the data for T = 30 
and 40 K are multiplied by factor 4 for clarity. From top to 
bottom, the data correspond to T =4.2, 10, 15, 20, 30 and 40 
K. The data measured on 36° bicrystals (9i w 36°, 62 ~ 0) 
at T =40 K in the same washer geometry are shown for com- 
parison (open circles). Right panel: The same for the sample 
No. 2. From top to bottom, the data correspond to T =4.2, 
10, 15 and 20 K. 

We assume that the period of a(Idc) at T =40 K, 
Aide =9.6 /iA, corresponds to A(p dc = 2tt. In order 
to determine the CPR we take ipdc = at a maximum 
or minimum of a. This is necessary in order to satisfy 
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I{ip = 0) = 0, as required by general principles a The ex- 
perimentally observed shift of the first local extreme of 
a(Idc) from Id c = (Fig. 2) can be due to flux trapped 
in the interferometer washer. Most probably, this flux re- 
sides in the long junction of the interferometer. The long 
junction does not play an active role because the Joseph- 
son penetration depth is much shorter than its length, 
and external fields producing I c i c are smaller than its crit- 
ical field. Nevertheless the long junction sets the phase 
difference for Id c = at the small junction. 
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Our experimental results can be understood as follows. 
Deviations from ideal geometry of the asymmetric 45° 
junction, (0\) = 45° + a\ and (62) = ct2, lead to a fi- 
nite value of I\. Thus, imperfections of the junction 
increase its critical current. For this reason we believe 
that samples Nos. 2-5 contain imperfections and from 
now we concentrate on nearly ideal junctions (such as 
sample No. 1) with |ai|,|a2| <C 1. For such junctions, 
the ratio I2/I1 exhibits the following temperature de- 
pendence. For T -► 0, \I 2 /h\ > 1. The region T ~ T c 
can be analyzed quite generally within Ginzburg-Landau 
theory. Let the electrodes be described by (macro- 
scopic) order parameters Ai^ = |A|e IV1,2 . Then the 
phase-dependent part of the energy of the junction is 
E = a[AiA£ + H.C.] + 6[(AhA£) 2 '+ H.C.] + . . . where 
a, b, . . . depend weakly on T.E-U Thus for T close to T c we 
estimate h cx |A| 2 cx {T c -T) and I 2 oc |A| 4 oc (T c -T) 2 , 
leading to \I2/Ii\ <C 1. These expectations are qualita- 
tively consistent with the experimental data shown in 
Fig. 4. 
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FIG. 3. a) Josephson current through the junction for the 
sample No. 1 as a function of the phase difference <p, deter- 
mined from the data in Fig. 2. The statistics of a(ip) was 
improved by folding the data back to the interval (0, n) and 
taking an average. The symmetry a(ip) = a(—tp) was as- 
sumed, b) The same for the sample No. 2. The data for the 
asymmetric 36° bicrystal at T — 40 K (open circles) is also 
shown. 

In Fig. 3, we show the CPR determined from the data 
in Fig. 2. For all curves we have performed a minimal 
necessary shift consistent with I(cp = 0) = 0. Thus, for 
the samples Nos. 1,2 we have assumed that at tp^c = 
a minimum of a((pdc) is realized. For an interferometer 
with a conventional s-wave weak link (and also for the 36° 
junction), at tpdc — a maximum of a(ipdc) is realized. 
Note that the maximum (minimum) of a(tpdc) at (fdc = 
implies a diamagnetic (paramagnetic) response of the in- 
terferometer in the limit of small applied fields. In Fig. 4 
we plot the coefficients I\ and I2 determined by Fourier 
analysis of the CPR for the sample No. 1 at various tem- 
peratures. With decreasing T, | X2 1 grows monotonically 
down to T = 4.2 K, while the I\ component exhibits only 
a weak temperature dependence. 



FIG. 4. Temperature dependence of the Fourier expansion 
coefficients Ji,2 determined from the experimental data in 
Fig. 3a. Solid lines are the Fourier expansion coefficients for 
the numerical data in Fig. 5. 

So far, our discussion was based solely on symmetry 
arguments. Let us attempt a more quantitative analy- 
sis of our data now. Two different microscopic pictures 
of asymmetric 45° Josephson junctions between d-wave 
superconductors have been considered in the literature. 
The first picture assumes a microscopically tetragonal 
material and an ideally flat interface.EirEJ Within this 
picture, there are two contributions to the Josephson 
current. The first is due to bulk states and in the tun- 
nel limit it iSj-well described by the Sigrist-Rice term 
I c in Eq. (||).0 The second is due to mid-gap states 
which develop close to the surfaces of unconventional 
superconductors.E£l I {(f) for the sample No. 1 calculated 
according to the model of Ref. [l2] is shown in Fig. 5. 
The experimental data can be fitted by a relatively broad 
range of barrier heights. However, if we require the 36° 
junction to be fitted by the same (or smaller) barrier 
height as for the 4SS junction, we conclude the barrier 
must be rather low .113 The T dependence of I{<f) requires 
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a choice of T c w 60 K in the non-selfconsistent theory 
of Ref. ||. The reduction from the bulk T c = 90 K is 
probably due to a combined effect of surface degradation 
and order-parameter suppression at the sample surface. 
The temperature dependence of the ratio of the tt and 
2-7T periodic components in I(<p) is seen to be in qualita- 
tive agreement with experimental data in Fig. 3a. This 
is explicitly demonstrated in Fig. 4 where we compare 
the experimentally obtained I\.2 with the results of the 
Fourier analysis of the curves in Fig. 5. The divergence of 
I2 as T — ► is an artifact of the ideal junction geometry 
assumed in Ref. [12[ If the finite roughness of the inter- 
face is taken into account, this divergence is cut off and 
the experimental data in Fig. 4 do indeed, resemble the- 
oretical predictions for a rough interfaces However, the 
nonselfconsistent theory of Ref. [l^ is unable to explain 
the experimentally observed steep CPR close to the min- 
ima of the junction energy (see Fig. 3a). In the limit of 
vanishing barrier height, the theoretical CPR does have 
steep portions, but these are located close to the maxima 
of the junction energy (see also Ref. O). 



Our present understanding of I(ip) in the asymmetric 
45° junction is only qualitative. There is considerable 
experimental evidenceffu that the grain boundary junc- 
tions are at most piecewise flat. However, we cannot say 
whether the shape of I(<p) is dominated by the mid-gap 
states in the microscopically flat regions, or by sponta- 
neous flux generation due to the spatial inhomogeneity 
of the junction. This issue requires further study. 

In conclusion, we have measured the magnetic field re- 
sponse of a single-junction interferometer based on asym- 
metric 45° grain-boundary junctions in YBa2Cu3C , 7_ a: 
thin films. Half-fluxon periodicity has been experimen- 
tally found, in agreement with theoretical predictions for 
d x 2_ y 2-wave superconductors. Hence, our results provide 
a novel source of evidence for the d-wave symmetry of the 
pairing state in the cuprates. 
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FIG. 5. I((p) calculated according to Eq. (64) of Ref. [L2| for 
a junction with 9i = 45.5°, 82 — 0, Ad = 1.5, k = 0.5, and 
T c = 60 K. /(</?) at T = 40 K for the 36° bicrystal (open 
circles) was calculated for the same parameters except for 
0i = 36°. 

In a different approach to the asymmetric 45° junction, 
one assumes a heavily twinned orthorhombic material 
(which is macroscopically tetragonal, hraKever) and/or 
a meandering interface with 0i = 0j(x)liJii3 Hence the 
critical current density j c (x) is a random function with 
a typical amplitude (\j c (x)\) ~ jc- If the average criti- 
cal current along the junction (j c ) <C j c , a spontaneous 
flux is generated in the junction, and \I%/Ii\ ^> lOEJ In 
particular, for (81) = 45° and {62) = 0, there is an equal 
amount of parts having positive and negative j c , leading 
t° (jc) — and I\ = 0. Note that also within the picture 
of Refs. [lj],[l8|, the d-wave symmetry of the pairing state 
is crucial, otherwise the condition (j c ) <C j c is difficult to 
satisfy. 
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